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Abstract 


Ozone measurements from ozonesondes, AROTAL, DIAL, and POAM III 
instruments during the SOLVE-2/VINTERSOL period are composited in a 
time- varying, flow-following quasi -conservative (PV-0) coordinate space; the 
resulting composites from each instrument are mapped onto the other instru- 
ments’ locations and times. The mapped data are then used to intercompare 
data from the different instruments. Overall, the four data sets are found to 
be in good agreement. AROTAL shows somewhat lower values below 16 km, 
and DIAL has a positive bias at the upper limits of its altitude range. These 
intercomparisons are consistent with those obtained from more conventional 
near-coincident profiles, where available. Although the PV-0 mapping tech- 
nique entails larger uncertainties of individual profile differences compared to 
direct near-coincident comparisons, the ability to include much larger numbers 
of comparisons can make this technique advantageous. 


Correspondence to: Leslie R. Lait (lait@code916.gsfc.nasa.gov) 
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1 Introduction 

The Kiruna deployments of the SOLVE-2 and VINTERSOL field experiments took place in 
January and February 2003. During this period, a number of different instruments measured 
stratospheric ozone. On board the NASA DC-8 aircraft were the Airborne Raman Ozone, 
Temperature, and Aerosol Lidar (AROTAL) and the Differential Absorption Lidar (DIAL), 
as well as in situ instruments such as FASTOZ and the Gas and Aerosol Measurement Sys- 
tem/Langley Airborne Measurement Spectrometer (GAMES/LAABS). Other instruments, 
such as the balloon-borne MkTV interferometer, flew on other platforms or took ground- 
based measurements. A number of sites launched ozonesondes in coordination with the 
campaigns. In addition, data from the Polar Ozone and Aerosol Measurement III (POAM 
EJ) solar extinction instrument on the SPOT-4 spacecraft were made available to the mission 
te ams 

Flying on board the NASA DC-8 together; the AROTAL and DIAL instruments were 
coincident, and the DC-8 flew over Ny Alesund in coordination with a number of the sonde 
launches there. But because many of the instruments operated on different platforms at dif- 
ferent times and locations, opportunities for measurement intercomparison were less than 
plentiful for most instruments. Aside from the AROTAL/DIAL data, the small number of 
near-coincident measurement sets make a statistical evaluation of inter-instrument differ- 
ences problematic. 

Other approaches are possible that can compare non-coincident data. The trajectory- 
mapping approach of Morris et al. (2000), the trajectory-hunting method of D anilin et al. 
(2003), and the MATCH technique of Rex et al. (1999) are examples. 

In this work, a quasi-conservative coordinate method is employed. Described in Schoe- 
berl and Lait (1991), this technique depends upon the premise that a reasonably long-lived 
trace gas should be well-mixed along contours of potential vorticity (PV) on a surface of 
constant potential temperature (i 9 ) (Leovy et al. , 1985). By using PV and 9 as coordinates, 
averaging mixing ratios near a set of points in that coordinate space should yield an accurate 
picture of a time-invariant trace gas distribution in PV-# space, in the absence of diabatic 
effects and chemical changes. For short time periods, these effects often may safely be 
ignored; for periods longer than 10 days or so, they must be taken into account somehow. 

PV-# analysis was used in Schoeberl et al. (1989) and Lait etal. (1990) to map measure- 
ments onto a three-dimensional field, and in Kyrd et al. (2000) and Lait et al. (2002) to 
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determine stratospheric ozone loss in the Arctic. The quasi-conservative coordinate method 
can be also useful in inter-instrument comparisons, as seen in Redaelli et al. (1994), Man- 
ney et al. (2001), and Randall et al. (2002). 

In this work, we intercompare ozone measurements from four instruments: AROTAL, 
DIAL, POAM III, and the ozonesondes. Section 2 describes each of these data sets briefly; 
then the analysis technique is described in section 3. Results follow in section 4. 

2 Data 

The SOLVE-2/V INTERSOL joint field experiment took place in January through early 
February 2003. We used the data from January 1 through February 10. 

The NASA Goddard Space Right Center’s AROTAL instrument is a lidar that uses 
Rayleigh scattering from Eximer and Yag lasers transmitting at 308 and 355 nm to mea- 
sure ozone, temperature, and aerosols. A more complete discussion of the instrument may 
be found in McGee et al. (2001) and Burris et al. (2002). For the SOLVE-2/VINTERSOL 
mission archive, AROTAL reports profiles every 22 seconds, averaged over 1.2 minutes. 
The altitude of the profiles depends on the altitude of the DC-8 aircraft, but over the middle 
of the flight it tends to range from approximately 14 km to 35 km. Vertical resolution of 
the reported data is approximately 150 m. Data were collected for 12 flights of the DC-8. 
To avoid problems with sunlight increasing noise in the measurements, only profiles taken 
where the local solar zenith angle is greater than 95° were used. 

A second lidar, the NASA Langley Research Center’s DIAL instrument, also flew on 
the DC-8. This instrument uses two YAG lasers transmitting at multiple frequencies to 
observe ozone below and above the aircraft. The DIAL data in the mission archive consist 
of profiles spaced about a minute apart. As with AROTAL, the altitudes covered change 
with the altitude of the aircraft, but typical coverage is from a few kilometers above the 
surface to around 25 km, with a small altitude gap near the aircraft itself. Vertical resolution 
reported is approximately 75 m. Data were collected for 14 flights (including two pre- 
mission test flights before the deployment to Kiruna). Details of this instrument may be 
found in Browell et al. (2003), Browell et al. (1998), and Richter etetal. (1997). 

A total of 213 balloon-launched sonde profiles from 21 ground stations were used in 
this analysis. These included special sondes launched for VINTERSOL, as well as those 
launched by the Meteorological Service of Canada, the World Meteorological Organization 
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network, Japan, and Russia. Data were used from the stations at Alert, Churchill, Eu- 
reka, Goosebay, Resolute, Stonyplain, Hohenpeissenbeig, Jokioinen, Sodankyla, Keflavik, 
Kuehlungsbom, Legionowo, Lerwick, Ny Alesund, Orland, Prague, Scoresbysund, Thule, 
Uccle, Salekhard, and Yakutsk. Altitude ranges vary greatly, but the sondes got as high as 
29 km. Reported vertical resolution also varies, from about 10 m to around 60 m. 

The PO AM III solar occultation instrument is described in Lucke et al. (1999). It is 
a nine-channel photometer that uses solar occulalion to measure atmospheric extinction in 
bands from 0.354 to 1.018 /mi to retrieve temperature and multiple species, including ozone. 
Fourteen profiles are taken each day, spaced around a latitude circle that moves slowly in 
time. Vertical resolution is approximately 1 km (Lumpe et al , 2003). 

To apply the quasi-conservative coordinate analysis, values of potential vorticity and po- 
tential temperature must be obtained at each measurement location and time. For this work, 
these are obtained by interpolating three-dimensional gridded analyses from the Data As- 
similation Office (now the Global Modelling and Assimilation Office) of NASA’s Goddard 
Space Flight Center. These analyses were chose because of their relatively fine horizontal 
and temporal resolution, as well as their altitude range. The analyses used were chosen from 
the “ final look” product generated from the GEOS-4 system. GEOS-4 was the successor to 
the GEOS-1 system documented by Pfaendtner et al. (1995). These data grids extend from 
1000 to 0.2 hPa, have a horizontal resolution of 1.25° longitude by 1.0° latitude, and are 
produced four times daily. 

3 Analysis 

In order to derive meaningful statistics for both coincident and non-coincident compar- 
isons, it was necessary to obtain uncertainties associated with the ozone measurements. To 
estimate these uncertainties for AROTAL and DIAL, a standard variance profile for each in- 
strument was constructed from the small-scale horizontal variations of the all profiles from 
all flights. For ozonesondes, we used the larger of 5% of the measured value, or the variance 
about a linear fit within a 5 km segment. Uncertainties in the POAM HI measurements were 
supplied in the data archive. 

The measurements, with their uncertainties, locations, and times were collected for each 
instrument (Note that for the sondes, the winds were used to estimate the horizontal loca- 
tion of each balloon during its ascent.) Modified potential vorticity (see Lait (1994)) and 
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potential temperature were interpolated from the meteorological analysis onto the measure- 
ment locations and times. 

By computing autocorrelations of the lidar data along the flight, we were able to estimate 
the minimum horizontal separation to ensure independence of profiles as approximately 
400 km for AROTAL measurements and 375 km for DIAL. These distances are roughly 
consistent with those used by Schoeberl et al. (2002). Because the profile sites for the 
sondes and POAM III data were widely separated, horizontal separations were not an issue 
for those instruments. 

Minimum vertical separations between a profile’s measurements were similarly obtained. 
For AROTAL, the vertical separation was estimated to be 3 km; for DIAL, it was 4 km. For 
the sondes, it was 3 km, and for POAM data, 4 km. 

Measurements closer than these horizontal and vertical separations were omitted from 
the analysis, so that the data going into the analysis could be considered independent of 
each other. 

The analysis itself is similar to that in Lait et al. (2002). A regular grid in a PV-0 coordi- 
nate space was constructed, and PV and 6 values interpolated from the analyses were used 
to locate each ozone measurement in the coordinate space. Data from January 1 through 
February 10 were used; hence, diabatic effects and chemical changes needed to be ac- 
counted for. These effects both show up as a change in ozone over time at a given point in 
PV-0 space. To first order, they can be dealt with by applying a weighted linear time fit to 
the data near a given P V-0 gridpoint. Each point was weighted inversely to its uncertainty 
and its distance from the P V-6 gridpoint being examined. 

This procedure yields a set of slopes and intercepts, one for each gridpoint in the PV-0 
coordinate space. From these, we can construct a composite field in PV-# space for any 
moment in time; moreover, given the meteorological analyses we can map this ozone field 
back into longitude-latitude-altitude space,. 

Note that for instrument comparison purposes here, only the overall evolution of the 
ozone field is relevant. Whether changes in that field are caused by diabatic effects or by 
chemical loss/production is of no concern, and so there is no need to try to separate those 
two effects. 

Data from two instruments are compared by constructing a PV-0 composite field from 
one instrument and mapping its ozone values onto the locations and times of the second 
instrument. The differences between the two ozone values and the uncertainties associated 


5 



with those differences are collected, and mean profiles of the biases are computed for the 
mission period taking the uncertainties into account 

4 Results 

To validate the analysis technique, several tests were applied. 

First, each instrument was compared ag ains t itself. That is, differences were character- 
ized between the original measurements and those from PV-# mapping of the same instru- 
ment’s data. This test should reveal any biases or distortions introduced by the analysis 
technique itself, and it should also reveal the degree to which noise is introduced by er- 
rors in the meteorological fields, departures from the well-mixedness assumption, and so 
on. Figure 1 shows an example of the self-comparison for AROTAL data The maximum 
difference, 0.1 ppmv near 22 km, is not statisticaly different from zero, and the rest of the 
profile is very close to zero. The self-comparisons for the other instruments show si milar 
results: very small average differences, with at most minor statistically insignificant fluctu- 
ations. 

The next test was to compute inter-instrument differences between AROTAL and DIAL 
using both near-coincident and non-coincident methods. Because both these instruments 
flew aboard the same aircraft, a large number of near-coincident profiles could be collected. 
For the 12 SOLVE-2 flights of the DC-8, the two closest profiles of the two instruments 
were chosen; each profile pair also had to be separated from all other profile pairs by at 
least 400 km. For each profile pair, the DIAL data were then interpolated to the AROTAL 
altitudes, and the two profiles were differenced. 

Figure 2 shows the differences and their average profile. Above 20 km, DIAL ozone 
values are systematically higher than AROTAL. Below 16 km, AROTAL values are higher. 

Comparing the AROTAL and DIAL data using the non-coincident PV-# analysis yields 
similar results (Figure 3). Note that the measured— reconstructed differences are consistent 
with the reconstructed— measured differences. Of course, because the data being compared 
were in fact coincident, this is no more rigorous a test of the noncoincident technique than 
the self-comparisons, evaluating the success of the next test. 

A more demanding test is to compare true non-coincident DIAL and AROTAL data. To 
accomplish this, AROTAL data from the even-numbered flights were compared with DIAL 
data from the odd-numbered flights, and vice-versa. The results, shown in Figures 4 and 5, 
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are consistent with the near-coincident comparisons, albeit with larger uncertainties (since 
they involved only half the data). Figure 4 matches the full data comparison quite well, 
while Figure 5 has greater uncertainties but is still roughly consistent with the others. 

Having confirmed that each instrument’s data compare well with themselves, and that the 
AROTAL-DIAL non-coincident comparisons are similar to the coincident comparisons, we 
proceeded to compare the other instruments’ data. 

Figure 6 shows the differences between AROTAL and the ozonesondes. There appears 
to be a slight bias below 15 km, with AROTAL being perhaps 0.3 to 0.4 ppmv higher than 
the sondes near 12-13 km. Above 25 km, the ozonesondes are fewer in number and their 
uncertainties are often larger, so that the error bars in the differences are much larger at 
those altitudes. 

Likewise, the differences between DIAL and the sondes is shown in Figure 7. Here, 
DIAL matches the sondes well at the lower altitudes, but DIAL is higher at the uppermost 
reaches of the instrument, above 25 km. Note also that the AROTAL-sonde and DIAL- 
sonde differences are consistent with the AROTAL-DIAL differences. 

Figures 8 and 9 show the POAM-AROTAL and POAM-DIAL differences, respectively. 
The POAM-AROTAL are qualitatively similar to biases noted in Lumpe et al (2003), where 
coincident comparisons were made between POAM III and AROTEL (an earlier version of 
AROTAL) and DIAL during the first SOLVE campaign in the winter of 1999-2000. The 
POAM-DIAL differences at the uppermost DIAL altitudes, however, are of different sign 
from those in Lumpe et al (2003) 

Figure 10 compares the ozonesonde data with the POAM III profiles. A small bias ap- 
pears near 20 km. At the uppermost altitudes, where the sonde measurements are fewer and 
less certain, there is also a suggestion of a small negative bias. 

These intercomparisons were repeated with meteorological analysis products from other 
institutions, and the results were similar. 

5 Conclusions 

Ozone measurements taken during the SOLVE- 2/VINTERSOL field experiment from four 
instruments (AROTAL, DIAL, POAM HI, and sondes) were compared. A quasi-conservative 
coordinate approach was employed to use non-coincident data for instrument interompar- 
isons. Several tests of the method were applied. First, each instrument’s data were self- 
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compared; the differences were zero, within the uncertainty associated with the technique. 
Second, the AROTAL and DIAL data were compared, with results similar to those from 
the straightforward near-coincident comparison. These tests demonstrate that the quasi- 
conservative technique does not introduce unusual or misleading artifacts into the data. 
To provide a true non-coincident data comparison, the DIAL data from the odd-numbered 
DC-8 flights were compared with AROTAL data from the even-numbered flights, and vice 
versa. These results were consistent with the other DIAL- AROTAL comparisons, although 
the uncertainties were larger. 

Finally, all four data sets were intercompared. The AROTAL ozone values were found to 
be 0.1 to 0.4 ppmv higher than DIAL and the sondes below 16 km, and DIAL values were 
0.3 to 0.8 ppmv higher than AROTAL and the sondes above 22 km. The ozonesonde data 
match well with POAM HI, except for a possible small bias in a region near 22 km. 

This PV-^analysis produces results with substantial uncertainties. The uncertainties as- 
sociated with the analyzed PV and 6 meteorological fields, limi ted sampling over regions 
of PV-f? space, potential failures in the assumptions necessary for the method’s validity 
(e.g., homogeneity around a circumpolar PV-0 tube), and even a breakdown in PV-ozone 
correlations at high altitudes and low latitude, can all contribute to the enlarged error bars. 
Certainly, then, direct comparison of large numbers of near-coincident measurements is 
preferred where it is possible. Nevertheless, the higher numbers of comparisons which are 
made possible by relaxing the requirement for near-coincidence, can improve the statistics 
so that the results are useful despite their uncertainties. 
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Figure Captions 


Fig. 1. Self-comparison biases for AROTAL data, expressed the difference between the original 
measurements and the measurements reconstructed through PV-0 mapping. Dots are die individual 
differences, the thin line is the mean difference profile, and the thick horizontal lines are the 90% 
confidence limits of the mean differences. The average measured ozone profile is shown in the right 
panel. 


Fig. 2. Near-coincident differences between AROTAL and DIAL. Gray dots repsesent individual 
differences, the solid line shows the mean difference profiles, horizontal lines show the 90% confi- 
dence limi ts of the mean profile, and the dotted lines show the standard deviation of the differences. 

Fig. 3. Differences between AROTAL and DIAL, using PV-0 anlysis. As in Fig. 1, the left panel 
shows the differences, and the right panel shows the average ozone profile. The mean differences 
(with their 90% confidence limits) are shown for the comparisons using AROTAL and DIAL as the 
data which is mapped onto the other instrument’s locations. 

Fig. 4. As in Fig. 3, except that AROTAL data from the even-numbered DC-8 flights are compared 
with DIAL data from the odd-numbered flights. 


Fig. 5. As in Fig. 3, except that AROTAL data from the odd-numbered DC-8 flights are compared 
with DIAL data from the even-numbered flights. 


Fig. 6. As in Fig. 3, except that AROTAL data are compared with data from the ozonesondes. 


Fig. 7. As in Fig. 3, except that DIAL data are compared with data from the ozonesondes. 


Fig. 8. As in Fig. 3, except that AROTAL data are compared with data from POAM HI. 


Fig. 9. As in Fig. 3, except that DIAL data are compared with data from POAM HI. 


Fig. 10. As in Fig. 3, except that ozonesonde data are compared with data from POAM III. 
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Figures 


Meas. AROTAL- Reconst. AROTAL Meas. AROTAL 



Fig. 1 . Self-comparison biases for AROTAL data, expressed the difference between the original 
measurements and the measurements reconstructed through PV-0 mapping. Dots are the individual 
differences, the thin line is the mean difference profile, and the thick horizontal lines are the 90% 
confidence limits of the mean differences. The average measured ozone profile is shown in the right 
panel. 


DC8 AROTAL- DC8 DIAL 



Ozone diff ppmv 


Fig. 2. Near-coincident differences between AROTAL and DIAL. Gray dots repsesent individual 
differences, the solid line shows the mean difference profiles, horizontal lines show the 90% confi- 
dence limits of the mean profile, and the dotted lines show the standard deviation of the differences. 
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Reconst AROTAL- Meas. DIAL 
Meas. AROTAL- Reconst DIAL 



Fig. 3. Differences between AROTAL and DIAL, using PV-# anlysis. As in Fig. 1, the left panel 
shows the differences, and the right panel shows the average ozone profile. The mean differences 
(with their 90% confidence limits) are shown for the comparisons using AROTAL and DIAL as the 
data which is mapped onto die other instrument’s locations. 


Reconst AROTAL (Odd)- Meas. DIAL (Even) 
Meas. AROTAL (Odd)- Reoonst DIAL (Even) 



Ozone diff (ppmv) 


Fig. 4. As in Fig. 3, except that AROTAL data from the even-numbered DC-8 flights are compared 
with DIAL data from the odd-numbered flights. 
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Reconst. AROTAL (Even)- Meas. DIAL (Odd) 
Meas. AROTAL (Even)- Reconst. DIAL (Odd) 



Fig. 5. As in Fig. 3, except that AROTAL data from the odd-numbered DC-8 flightsare compared 
with DIAL data from the even-numbered flights. 


Reconst. AROTAL- Meas. Sondes 
Meas. AROTAL- Reconst. Sondes 



Ozone diff (ppmv) 


Fig. 6. As in Fig. 3, except that AROTAL data are compared with data from the ozonesondes. 
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Reconst POAM3- Meas. DIAL 
Meas. POAM3* Reconst DIAL 



Fig. 9. As in Fig. 3, except that DIAL data are compared with data from POAM HL 


Reconst PQAM3- Meas. Sondes 
Meas. PQAM3- Re co n st Sondes 
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Ozone diff (ppmv) 


Fig. 10. As in Fig. 3, except that ozonesonde data are compared with data from POAM HL 
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Reconst. DIAL- Meas. Sondes 
Meas. DIAL- Reconst. Sondes 



V 

Fig. 7. As in Fig. 3, except that DIAL data are compared with data from the ozonesondes. 


Reconst. POAM3- Meas. AROTAL 
Meas. POAM3- Reconst. AROTAL 



Fig. 8. As in Fig. 3, except that AROTAL data are compared with data from POAM IQ. 
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Popular Summary of 

“Non-coincident Inter-instrument Comparisons of Ozone 
Measurements Using Quasi-conservative Coordinates” 
by Lait, et ai. 


This paper compares stratospheric ozone measurements from four instruments during 
the SOLVE-2/VINTERSOL field experiment in January-February 2003. The four instru- 
ments are: the NASA Goddard Space Flight Center AROTAL airborne lidax, the NASA 
Langley Research Center DIAL airborne lidar, the satellite-based POAM III instrument of 
the U. S. Naval Research Lab, and various ozonesondes launched during the period. Be- 
cause the number of coincident measurements was small for most instruments, a special 
non-coincident technique was employed here that uses conserved meteorological quantities 
to map data from one instrument onto another instrument’s measurement locations. We find 
that most instruments agree reasonably well over most of their altitude ranges. AROTAL 
has a somewhat positive bias below around 13 km, and DIAL has a stronger positive bias 
above 22 km. 
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